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Abstract 

Aims. While chemical composition spreads are now believed to be a universal characteristic of globular clusters (GCs), not all of them 
present multiple populations in their color-magnitude diagrams (CMDs). Here we present a new scenario for the formation of GCs, 
in an attempt to qualitatively explain this otherwise intriguing observational framework. 

Methods. Our scenario divides GCs into three groups, depending on the initial mass (Mi) of the progenitor structure (PS), as follows, 
i) Massive PSs can retain the gas ejected by massive stars, including the ejecta of core-collapse SNe. ii) Intermediate-mass PSs can 
retain at least a fraction of the fast winds of massive stars, but none of the core-collapse SNe ejecta. iii) Low-mass PSs can only retain 
the slow winds of intermediate-mass stars. Members of the first group would include a> Centauri (NGC 5139), M54 (NGC 6715), M22 
(NGC 6656), and Terzan 5, whereas NGC 2808 (and possibly NGC 2419) would be members of the second group. The remaining 
GCs which only present a spread in light elements, such as O and Na, would be members of the third group. 

Results. According to our scenario, the different components in a) Cen should not display a sizeable spread in age. We argue that 
this is consistent with the available observations. We give other simple arguments in favor of our scenario, which can be described in 
terms of two main analytical relations: i) Between the actual observed ratio between first and second generation stars (R^q) and the 
fraction of first generation stars that have been lost by the GC (S£); and ii) Between S L and Mj. We also suggest a series of future 
improvements and empirical tests that may help decide whether the proposed scenario properly describes the chemical evolution of 
GCs. 

Key words, globular clusters: general; globular clusters: individual: NGC 2808, NGC 6715, NGC 5139, NGC 6656, Terzan 5; stars: 
abundances; stars: formation; stars: evolution 
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1. Introduction 

For many years it was believed that globular clusters (GCs) fol- 
low with excellent precision the canonical zero-order approxi- 
mation for the formation of simple stellar populations, where a 
chemically homogeneous gas cloud collapses under the action 
of its own potential well forming stars of different masses at the 
same time. However, recent observations have shown that this 
hypothesis is quite far from reality. 

One of the peculiarities observed is a spread in the abun- 
dance of light elements in stars of all GCs studied to date. This 
phenome non is very well rep resented by the O-Na anticorrela- 
tion (see iGratton et al.l I2004 for a thorough review about the 
chemical composition of GC stars). This anticorrelation differs 
in detail from one GC to the next, and its extension seems to de- 
pend on the present-day GC mass, but it does not seem to depend 
on the cluster's metallicity (ICarretta et al. 2010a). 

Due to the fact that massive GCs have deeper potential 
wells, it is commonly thought that the spread of Fe-peak ele- 
ments would only be present in the most massive GCs. However, 
spectroscopic evidence has revealed heavy element variations in 
GCs which are less massive than others which do not show a 
spread in these elements. This is the case of M22 (NGC 6656) 
and NGC 1851, with masses of ~ 5.4 x 10 5 M o , and ~ 5.6 x 
10 5 M o , respectively, which show heavy element variations, 
while other GCs (such as the extensive studied M3, with a mass 



~ 8.4 x 10 5 M o ) do not show this peculiarity^ In fact, while 
M22 shows an iron spread between [ Fe/H| » -1.9 and —1.6 , 
as well as sp reads in other element s (iPilachowski et all 1 1982: 
[ Lehnert et ail I199U iDa Costa et ail 120091: iMarino et all [2009, 
1201 lbl) . NGC 1851 harbors two groups of stars with different 
[Ba/Fe] ratios, but, depending on the authors, with a uniform 
[Fe/H] (lYong & Gmndah1ll200a IVillanova et a U l20Toh or with 
a spread in metallicity of A[Fe/H] ~ 0.2 dex ( Carretta et al.l 
l2010bll201 ll who also associate NGC 1851 to a merger of two 
separate GCs). From a deep study of their color-magnitude dia- 
grams (CMDs), these GCs show that the simple stellar popula- 
tion hypothesis is r uled out, as revealed by an observed sub giant 
branch (SGB) split dMarino et alj|2009l iMilone et alj|2008l) . 

Another GC presenting an SGB split is NGC 6388 dPiottol 
2009), with a mass of 1.4 x 10 6 M.. ICassisi et alJ (I2008T) have 
pointed out that such SGB splits could be related to an en- 
hancement in CNO abun dances, or an in t ernal age difference 
of a few Gyr. However, Villan ova et alJ d2010h have shown 
that for NGC 1851 the C+N+O content would be constant 
between the two groups with different [Ba/Fe]. Moreover, 



Di Criscienzo et al. (2010) have pointed out that such SGB splits 
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1 Here, as in the rest of this pa per, we estimate the mass of each GC 
using M v from the Harris ( 1996, revision 2010) catalog and a mass-to- 
light rat io M/Ly = 3 MpL: 1 , as p redicted by simple stellar population 
models dBruzual & Chariot 2003). 
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may also be caused by variations in the helium abundance, 
though only at moderately high metallicity. 

For the moment, NGC 2808 (1.4 x 1O 6 M ) seems to be 
a fairly special case, with at most a small spread in metals 
( TFe/Hl = —1.10 ± 0.03) but a wide O-Na anticorrelation 
dCarretta et al.ll2006b . Indeed, this GC has three populations ob- 
served down to the main sequence (MS), which suggests differ- 
ent helium abundances, ranging from Y _ 0.25 to Y = 0.3 7 
dNorrisI 120041 ID' Antona etaD 120051: IPiotto et alJ 120051 120071) 
This l ast peculiarity was in fact predicted by ID' Antona et al.l 
(120021) from a study of the cluster's horizontal branc h (FIB 
morphology, and recently tested by Dalessan dro et al.l d2011l 
using high-resolution far-UV and optical images, as well as 
Pasquini e t al.l (1201 ll) using visual and near-IR spectroscopy. 
Another GC which may share properties similar to NGC 2808 's 
is NGC 2419, wit h a present-day mass o f 1.4 x 1O 6 M 
dDi Criscienzo et alJ 1201 it ICaloi & D'Antonal 120111: see also 
ICohen et alj|2010l) . 

From photometry, the case of a> Cen, the most massive 
Galactic GC (3 x \0 b M n \ has been known for a long time dGeverl 
ll967tlCannon & Stobid 19731) . although only recent studies have 
unveiled some of its most intriguing properties. This GC has 
at least three separate, well-defined red giant branches (RGBs), 
an extended HB morpholog y, and a large number of subiumi- 
nous e xtreme HB stars (e.g . . iD'Cruz et al.lll996t IWhitnev et al. | 
1998; iPulone et al.l 1 19981: IrrCruz et al.l l2000t iPancino et al. 
2003 : iRev et al.l l2004t iFrevhammer et al l 120051: ICassisi et al.l 
2009; iBellini et al.l 120091: ICalamida et al.l I2009L (201 Oh . Recent 



Hubble Spac e Telescope (HST) observations have also shown 
a triple MS dBedin et all 12004 iBellini et al.l 120091 1201 Oh . in 
addition to at least four SGBs ( Vill anova et alJ 120071) . From 
spectroscopic studies, this GC also shows a complex behav- 
ior. Based o n spectroscopically derived chemical compositions 
of MS stars, IPiotto et aTl (120051) have suggested that the bluest 
MS is highly enriched in helium (Y = 0.35), while metallic - 
ities of SGB and RGB sta rs show a wide distribution, from 
rFe/Hl ~ -2 to ~ -0.4 dHilker et al.l 12004 I Villanova et al.l 
120071: Ijohnson et al.l 12009: Johns on & Pilachowskil 1201 d and 
references therein). These different metallicity groups present 
some intriguing properties, including a smaller spread in the 
abundances of some elements (e.g., C, N, Ca, Ti, and Ba; see 
Fig. 16 in Villanova et al. 2007) at the more metal -rich end, 
compared to metal-poor stars in the cluster. It would thus ap- 
pear that metal-rich stars were formed all at the same time, from 
a fairly chemically homogeneous cloud. Metal-poor stars may 
also have formed all at the same time, but from a chemically 
/nhomogeneous cloud. Alternatively, these metal-poor stars may 
have formed at diffe rent times, with m aterial processed in differ- 
ent ways (see also lGratton et al.ll201 lh . 

Another GC with properties similar to a> Cen's, albeit at a 
much less extreme level, is M54 (2.3 x 10 6 M H ), in the nu cleus 
of the Sagittarius dwarf galaxy (e.g.. ISiegel et al.l 120071) . The 
fact that M54 is the second most massive GC in the Galactic 
neighborhood strongly suggests that u Cen too may be the rem- 
nant of a formerly much larger structure (e.g., lDinescul f2002: 
lAltmann et ail 120051: IWvlie-de Boer et al.ll201o[ and references 
therein). Other GCs where internal dispersions in the Fe-peak 
elements have recently been detected, and which may thus 
share some similarities in th eir chemical evolution with u> Cen 
and M54, inc lude M22 (e.g. jDa Costa et al.ll2009l iMarino et al 



2009.1 201 lbh and Terzan 5 dFerraro et al.ll2009t lD' Anton a et al 
2010tlOrigliaetal.ll2011l) . 

An intriguing GC is 47 Tucanae (NGC 104), which in spite 
of being one of the most massive Galactic GCs (1.4 x 1O 6 M ) 



presents but a moderately wi d e SGB - a feature that is inter- 
preted by iDi Criscienzo et aD (1201 Ol) as being due to the pres- 
ence of subpopulations with slightly different helium abun- 
dances, instead of either an internal spread in age or C+N+O. 
Importantly, the lack of prominent subpopulations in 47 Tuc sug- 
gests that the present-day GC mass does not correlate tightly 
with the presence and extent of the multiple populations ob- 
served in GC CMDs, like ly indicating important mass evolu- 
tion in these objects (e.g., Vesperinilll998tlFall & Zhangfe OQ 1 ; 
Bekki & Nomsl2006LlLV Antona et alJ2007tlMcLaughlin & Falll 



20081: Conrovl 20lH "and references therein). On the other hand, 
Nata f et alJ ( 201 1 ) have very recently found a radial gradient in 
the luminosity of the RGB bump in 47 Tuc, interpreting this as 
evidence of a second generation of stars forming deep within the 
cluster's potential well. 

In this paper we present a scenario which tries to explain 
qualitatively the observed CMD peculiarities, where the main 
parameter that defines the present-day CMD is the initial mass 
(M/) of the progenitor structure (PS), as opposed to the GCs 
present-day mass (Mcc)- We assume that all GCs have lost an 
unknown amount of mass which is not necessarily related with 
Mqq, but which can be estimated assuming a simple analytical 
model. 

Our paper is structured as follows. First, we review in §|2] 
some of the recent scenarios for the formation of GCs. This 
is followed in §|3]by an explanation of our suggested scenario. 
Then in §|4] we estimate the M/ value of the PS which in our 
scenario has given birth to the present-day GC, based on the ob- 
served ratio between second and first generation stars (R s ^)- In 
this section we also discuss some open problems with our sce- 
nario. Finally, in $5] we summarize our main results. 

2. Previous Scenarios 

For decades there have been many attempts to establish a 
scenario to explain the formation and ev o lution of GCs (e.g., 
iMatsunami et al.l 119591: iPikel'Nerl 119761: IPiFazio & Re nzini 



1980; Pi Fazi ol Il986t henorio-Tagle et alJ Il986t iBrown etal 



1994 



1991; Murr ay & Lin[ 
Richtler & Fichtnerl 119931: 
19951: ISalaris & Weissl 1 19971 



IPrice & Podsiadlowskil I ! 993 
Murra y & Linl 119931: lBrown ~ 
iNakasato et al.l 120001) . but only 



[993b 
et al.1 



recently has it become clear that the common assumption that 
they are simple stellar populations (SSP) is generally incorrect - 
and, importantly, not only in the well-known cases of to Cen 
and (more recently) M54. In fact, there is increasing evidence 
that virtually all GCs show some level of chem ical composition 
inhomogeneity (e.g., ICarretta et al.l l2009albh . which in some 
cases may lead to multiple sequences in the observed CMDs. 
In what follows we summarize some of the scenarios that have 
rece ntly been prop osed to explain the new empirical evidence 
(see [Renzini 200J, for a more detailed comparison). 



2. 1. \Decressin eTa/1 \2Q07a) Scenario 

iDecressin et al.l d2007al) have presented a scenario to explain 
the multiple populations observed in some massive GCs which 
invokes pollution fro m fast rotating massive stars (FRMS, 
IDecressin et al.ll2007bl) . The latter are stars that have reached the 
critical rotation velocity, and have formed an equatorial gas disk 
around them. Unlike the normal winds of massive stars, which 
have velocities ranging from hundreds to a few thousand km/s, 
the equatorial disk loses mass with velocities lower than 50 km/s, 
which can then easily be retained by the potential wells of GCs. 
In this scenario, FRMS ejecta are mixed with pristine gas to 
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form second generation (SG) stars near their massive progeni- 
tors. Even though this scenario may help explain the O-Na an- 
ticorrelation observed in all GCs, the multiple MSs observed in 
some GCs cannot be explained without invoking discrete helium 
abundances, whereas this scenario might more naturally lead t o 
a smooth spread in Y, as was also pointed out bv lRenzinl ((2008). 

2.2. \D'Ercole ef a/1 \2008) Scenario 

In Id] Ercole et"al1 d2008l) . the formation of multiple populations 
is attributed to the g as ejected by intermediate-mass AGB (or 
"super-AGB") stars (iD'Antona et alJfl983l: IVentura et al.ll200ll 
2002; iThoul et al.l l2002), with this gas starting to collect in the 
GC core after the last core-collapse SN explosion. In this sense, 
the formation of a cooling flow is studied by the authors using 
one-dimensional hydrodynamical simulations. Such a cooling 
flow collects the gas expelled by super-AGB stars that is used to 
form the SG stars. Using N-body simulations, the authors study 
the dynamical evolution of both populations, showing that if SG 
stars are formed in the center of the GC, most of the expelled 
stars will be first generation (FG) stars, which can result in a GC 
dominated by SG stars. 

For their "standard" scenario, hydrodynamical simulations 
show that the star formation rate increases since 10 Myr (the 
time of the last core-collapse SN) until 40 Myr (the time of the 
first SN la). However, the effect of core-collapse SNe belonging 
to the SG, which must stop the formation of SG stars after a few 
Myr, is not considered. 

Note also that every SG star that does not explode as a core- 
collapse SN is a possible progenitor of a type la SN. In fact, the 
more massive the progenitor, the closest the remm aining white 
dwarf will be to the C handrasekhar mass (e.g., ISalaris et al.l 
l2009HKaliraTet al. 2009). Therefore, after the core-collapse SNe 
period, there must also be a period of type la SNe, which will 
inhibit the formation of stars with gas ejected by more massive 
intermediate-mass stars. Naturally, the frequency of these events 
must decrease with time, instead of remaining constant. In our 
opinion, this poses considerable problems for the explanation of 
multiple populations with high helium abundances in GCs using 
intermediate-mass stars. 

2.3. \Conrov & Sperael \2011\) Scenario 

In the IConrov & SpergeJ d201 ll) scenario, the FG stars are 
born in an intracluster medium previously enriched in metals 
(present-day GC metallicity), while after some Myr the remain- 
ing gas is completely expelled by core-collapse SNe. The GC is 
then ready to start the accretion of the intermediate-mass AGB 
ejecta for several xlO 8 yr. During this time, the GC is also ac- 
creting pristine material from the ambient interstellar medium 
(ISM), which is mixed with the AGB ejecta to form SG stars. 
Finally, SG core-collapse SNe, and later FG type la SNe, begin 
to explode, keeping the GC gas-free, which stops the star forma- 
tion process. 

Even though the accretion of material from the ISM is plau- 
sible, it is unlikely that the chemical composition of this ac- 
creted material is similar to the GC itself. Moreover, if core- 
collapse SN explosions cleaned the intracluster medium, it must 
also clean the surrounding ISM or/and increase its metallicity, 
implying that SG stars will be formed with different metal abun- 
dances. 

To close, we note that iD'Ercole et alj d201 lb have very re- 
cently also argued against this scenario, pointing out in particular 



that it is not possibl e to form very O-poor s t ars through the pro- 
cess envisaged bv lConrov & Spergeld201 lh . lD'Ercole et alJalso 
point out, following iRenzinil d2008l) . that the Conrov & Spergel 
scenario also fails to account for the high amount of He that is 
needed to explain the highly He-enriched populations that ap- 
pear to be present in some GCs (see ©. 

2.4. \Carretta et al\ \2010a\) Scenario 

ICarretta et al.l d2009bl and references therein) have made a re- 
markable effort to obtain a homogeneous spectroscopic database 
for 17 GCs, where they have focused mainly on the O-Na anti- 
correlation. These data were used to divide the stellar popula- 
tions in GCs into three components: i) The primordial popula- 
tion, with stars with O and Na abundances similar to field stars 
(i.e., stars with [Na/Fe] < [Na/Fe] min + 0.3, where [Na/Fe] min 
is the lowest Na abundance detected in that GC); ii) The inter- 
mediate population, with stars with [O/Na] > -0.9 that do not 
belong to the primordial population, amounting to 60% to 70% 
of the observed stars; and iii) The extreme population, with stars 
with [O/Na] < -0.9, which are not presented in all GCs. 

Using these data and other global paramete rs characterizing 
the GCs in their sample, ICarretta et al.l (1201 Oal) have proposed a 
scenario where GCs were formed in three different stages. First, 
a precursor population of stars was formed when the unborn GC 
(with a size of ~ 100 pc, and made of gas and dark matter) in- 
teracts strongly with other structures. Core-collapse SNe of this 
population enrich the remaining gas, and trigger the formation of 
the primordial population. Then, the gas ejected by primordial 
FRMS or super-AGB stars give rise to a gas cloud chemically 
enriched in the center, where the second generation (SG) of stars 
is born. Finally, SG core-collapse SNe clean the remaining gas, 
thus halting star formation. During this time, the structure has 
lost all its dark matter content, almost all the precursor stars, and 
a large fraction of the primordial population. 

While this appears to provide a promising framework, there 
are some points which deserve a deeper inspection. Primordial 
low-mass stars, which are not observed in the actual GC be- 
cause they were expelled from the initial structure, must still 
be present in some place. Therefore, according to this sce- 
nario, a large number of metal-poor stars must be present in 
the field. However, to date there are only 174 known stars 
with fFe /Hl < -3, and 65 9 stars with [Fe/H] < -2 (SAGA 
database^ Suda et al. 2008). Alternatively, the initial mass func- 
tion (IMF) must be top-h eavy at low metallicities (Skillman 
2008; iKomiva et all l2010h . Note, in addition, that if the gas 
ejected by the precursor core-collapse SNe is retained in the ini- 
tial structure, the ejecta of the primordial core-collapse SNe must 
also be retained, which would increase the metallicity of the SG 
stars unless the initial structure had already lost a large fraction 
of its mass in a short period of time. 

2.5. A Final Objection Against the Super-AGB Scenario for 
the Population Richest in He 

D' E rcole et al. I d20Toh have shown that, in order to reproduce 
the observed O-Na anticorrelation, as well as the suggested He 
spread in NGC 2808, SG stars must form from super-AGB ejecta 
only, whereas in the case of M4 (NGC 6121) SG stars must be 
formed with a large amount of pristine gas mixed in. This is 
very surprising, since NGC 2808 is ten times more massive than 
M4 (1 .8 x 1O 5 M ), whereas it should be easier for a massive PS 

2 http://saga.sci.hokudai.ac.jp 
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Figure 1. Schematic representation of our scenario for low-mass PSs. Stars of the first and second generations are represented as 
red and orange, respectively. Sizes indicate the different star masses, where large, medium, and small circles represent massive, 
intermediate-mass, and low-mass stars, respectively. Arrows show the gas flux, with the arrow color indicating the chemical com- 
position following the same color scheme as for the stars of the different generations, and the arrow size being roughly proportional 
to the gas speed. Colored areas represent zones of gas accumulation. In each panel an approximate time for these events is given. 
Panels represent: a) Formation of FG stars; b) Pristine gas expulsion due to winds of FG massive stars; c) Core-collapse SN ejecta 
are not retained in low-mass PSs; d) Gas accumulation from winds of FG super-AGB and AGB stars; e) Formation of SG stars; 
f) SG core-collapse SNe period; g) Present-day O-Na anticorrelation. 



than for a low-mass PS to retain pristine gas. In this case, M4 
would have lost many more stars since its formation than did 
NGC 2808, implying a more massive PS in the past - but not 
so massive as to have retained any (unobserved) trace of met- 
als from SN explosions. However, due to their deeper potential 
wells, massive PSs are expected to lose a smaller number of stars 
than less massive PSs. This is contrary to what would have been 
expected if M4 were related with a massive PS but NGC 2808 
with a less massive one, unless the M4 progenitor has been af- 
fected by stronger dynamical interactions in the course of its life- 
time. 

Finally, we believe it is important to emphasize that, in our 
scenario, super-AGB and/or normal AGB stars do play an impor- 
tant role in the formation of multiple populations. In fact, as we 
shall see, while our scenario implies that these stars cannot be 
the progenitors of the most helium-rich population observed in 
NGC 2808 and u> Cen, it does also indicate that they are likely to 
be the progenitors of the stars with helium abundances falling in 
between the "normal" (i.e., non-He-enriched) and most He-rich 
stars. 



3. Our Scenario 

Our new scenario is divided in two stages, and is schematically 
described in Figs. [T] [3] and[5]l. The first stage gives the basis for 
the formation of any GC, reflecting the conditions in the ISM 
that allow a gravitationally bound "proto-GC" to form (panels 
labeled "a" in Figs.[TJ[3]|5]l. In turn, the second stage depends on 
the PS's initial mass M/, and can have three different outcomes. 
More specifically, i) a low-mass PS only retains the mass ejected 
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Figure 2. CMD of a low-mass PS after 13 Gyr. PGPUC 
isochrones represent the two surviving populations for this sys- 
tem: FG (red lines) and SG (orange lines). 

via slow winds from intermediate-mass stars; ii) an intermediate- 
mass PS does not retain the core-collapse SNe ejecta, but does 
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Figure 3. As in Fig. \T\ but for intermediate-mass PSs. Here the first, second, and third generation of stars are represented in red, 
blue, and green, respectively. Panels represent: a) Formation of FG stars; b) Gas accumulation from winds of FG massive stars and 
pristine gas; c) First core-collapse SN explosions, which trigger the formation of SG stars and expel the gas that has not arrived to the 
core; d) FG and SG SN explosions; e) Gas accumulation from winds of super-AGB and AGB stars of the FG and SG; f) Formation 
of TG stars; g) Present-day O-Na anticorrelation. 



retain at least a fraction of the wind from massive stars; and iii) a 
massive PS retains the winds of massive stars, and also the ejecta 
of core-collapse SNe. 



3.1. First Stage 

1. As usual, the first step for the formation of a GC begins 
with the gravitational collapse of a cloud with the present- 
day metallicity of the GC. 

2. FG stars are then formed following a homogeneous distribu- 
tion throughout the GC. They are initially embedded inside 
the ISM gas that was not used up to form stars (panels la- 
beled "a" in Figs.[T|[3] and|5]l. 

3. Assum ing a star forma tion efficiency (e) between 20%, and 
40% (Par mentier et a l. 2008), the remaining gas has a mass 
between 0.8 M/ and 0.6 M/, but is distributed over a larger 
volume than is required for triggering a new local fragmen- 
tation process to form a new generation of stars. This implies 
a decrease in the local gas pressure for the same potential 
well, with the consequence that the remaining gas starts to 
fall again into the center of the cluster. 

4. The final step of this stage begins when massive stars eject 
their envelopes at high velocities, which collide with the 
falling gas. This leads to a decrease in the speed of the out- 
bound massive star ejecta, and also to a decrease in the speed 
of the inbout gas (panels labeled "b" in Figs. Q] [3] and|5]l- 
The chemical composition of the massive star ejecta is ini- 
tially quite similar to the primordial one, however after a 
while (depending on the initi al stellar mass) this ejec ta con- 
sists of almost only helium (Limongi & Chieffi 2007), with- 
out differences in heavy element abundances except mainly 
for those elements that participate in the CNO, NeNa, and 
MgAl cycles. 



3.2. Second Stage for Low-Mass PSs 

For low-mass PSs, the second stage takes place as follows (see 
also panels a through f in Fig.Q]): 

5a. Low-mass PSs are unable to retain the FG massive star ejecta 
because of their shallower potential well, and the consequent 
low velocity of the infalling primordial gas. For the same 
reason, the PS center does not have enough material to trans- 
form all the kinetic energy of massive star ejecta into thermal 
energy. The formation of a viable star-forming cloud in the 
core is thus inhibited (panel b in Fig.[TJ. 

6a. The next step is driven by the FG core-collapse SNe explo- 
sions, which clean completely the PS of the remaining pri- 
mordial gas (panel c in Fig. |T}. 

7a. In this case, SG stars are formed only with (properly diluted) 
gas ejected by super-AGB and/or AGB stars (panels d and 
e), with the chemical composition of the ejecta depending in 
detail on the stellar mass. 

8a. This is potentially a continuous process, with renewed 
cleansing of the intracluster gas taking place after each new 
star formation event (panel f in Fig.Q]). 

If the above scenario is correct, each successive stellar gen- 
eration would leave a different signature in frequently observed 
diagnostic planes, such as the [Na/Fe] vs. [O/Fe] diagram. This 
diagram is schematically shown in panel g of Figure Q] The an- 
ticipated impact on the HR diagram is shown in Figure |2l where 
our new set o f Princeton-Go ddard-PUC (PGPUC) isochrones 
has been used dValcarce etaDl20TTh . 

3.3. Second Stage for Intermediate-Mass PSs (NGC 2808's 
Progenitor) 

A schematic representation of the second stage of the formation 
of intermediate-mass PSs, which are related to the progenitors of 
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Figure 4. As in Figure [2] but for an intermediate-mass PS. 
PGPUC isochrones represent the different populations presented 
for this PS: FG (red lines), SG (blue lines), and TG (green lines). 

GCs like NGC 2808 (and possibly NGC 2419), is shown in pan- 
els b through f of Figure[3] and can be summarized as follows: 

5b. Due to the fact that intermediate-mass PSs have deeper grav- 
itational potential wells than do low-mass PSs, the falling 
gas reaches a higher speed, and a fraction of the massive star 
ejecta is retained before their progenitors explode. However, 
in the outer part of the PS the massive star ejecta are try- 
ing to escape, while the pristine gas is infalling (panel b in 
Fig. [3}. Additionally, in the PS core, where the primordial 
gas has been accumulating, the gas ejected by massive stars 
is slowed down, transforming the kinetic energy into thermal 
energy, which delays further star formation. 

6b. When FG core-collapse SNe explode, their ejecta compress 
the center cloud, thus triggering the SG's star formation. 
These SG stars are highly enriched in helium due to the mas- 
sive star ejecta, but at the same time they are not heavily 
enriched in metals. Assuming that SN explosions are nearly 
symmetrical and that their precursors are not too close to 
the center, only a small fraction of the SN ejecta, which is 
metal-enriched, will be mixed with the core cloud. This event 
also completely removes the outer mixed gas from the clus- 
ter (panel c in Fig. |3j, since the gravitational potential well 
is not deep enough and the mass of the outer falling gas is 
insufficient to retain the SNe ejecta. Thus, TG stars will ac- 
cordingly not be metal-enriched. 

7b. Now, the ejecta of massive SG stars are not retained in the 
case of an intermediate-mass PS, because the infalling gas 
is only produced by intermediate-mass FG stars. If any such 
gas is initially retained, it will eventually be expelled by SG 
core-collapse SNe or FG type la SNe explosions, or will 
form but very few stars (panel d in Fig. [3). 

8b. After this second cleansing of the cluster, a new cloud is 
starting to form in the cluster center, using the intermediate- 
mass stellar ejecta (mass lost at low velocity) from FG and 
SG stars (panel e in Fig. [3). The chemical composition of this 



new cloud falls between both generations, as a consequence 
of the slope of the IMF (which favors low-mass stars) and 
the mass ratio between both generations. Here, the first stars 
of the third generation (TG) will be created. 
9b. Such star formation and cluster ISM cleansing stages con- 
tinue, with each successive stellar population being chemi- 
cally more similar to FG stars, while at the same time less 
numerous (panel f in Fig . [3} . 

In panel g of Figure [3] we show schematically the expected 
shape of the diagnostic [Na/Fe] vs. [O/Fe] diagram. The corre- 
sponding HR diagrams for this class of PS are shown in Figure|4] 

3.4. Second Stage for Massive PSs (<x> Cen's Progenitor) 

Figure [5] is a schematic representation of our scenario for mas- 
sive PSs, which is explained in the following paragraphs. Here, 
one can identify to Cen as a prototype of massive PSs - but other 
posssible examples include M22, M54, and Terzan 5 (see ©. 

5c. As in 5b, in the PS core a cloud has been formed that is 
highly enriched in helium. Its material comes from massive 
star ejecta and a fraction of the remaining gas that was not 
used to form FG stars. However, in the outer part of the GC 
the massive star ejecta is trying to escape, while the pris- 
tine gas is infalling with a velocity (and a total mass) that is 
higher than in the case of intermediate-mass PSs. 

6c. The first core-collapse SNe explosions compress the core 
gas, which triggers the SG star formation episode (panel c 
in Fig. [5]). As in 6b, these SG stars are highly enriched in 
helium (without however having significantly higher metal 
abundances). However, contrary to what happens in the case 
of intermediate-mass PSs, the higher potential well does al- 
low the gas in the outer part of the cloud to be retained. Most 
of the SN ejecta tries to escape the cluster, merging with the 
falling gas in the process. This event efficiently mixes both 
gas components, and delays the moment of arrival of the fur- 
ther mixed gas to the core. 

7c. After a while, the highly metal-enriched material is mixed 
with the gas which was not used to form SG stars, creating 
a new cloud in the core of the cluster (panel d in Fig. [5j. 
This cloud is also fed by SG massive stars and by both mas- 
sive and intermediate-mass FG stars - and these provide the 
chemical ingredients that will characterize the cluster's third 
generation (TG) of stars (panel e in Fig. [5). Unfortunately, 
this population has a helium abundance which is difficult to 
predict, because of the several sources of material that are in- 
volved. However, we expect that TG stars have some degree 
of helium enrichment. 

8c. Now, with three stellar generations in the cluster, the pro- 
cess of star formation continues, but each time with ma- 
terial mainly processed by less massive stars (super-AGB 
and/or AGB stars). As a result, the newborn stars belonging 
to the fourth generation (4G) will have a chemical compo- 
sition which is a mix between the three preceding genera- 
tions. However, in our scenario each new star that is formed 
will have a chemical composition more similar to that of 
the primordial gas, owing to two main factors: i) The IMF 
slope, which implies more stars with low masses than mas- 
sive stars, and ii) The relatively small amount of mass used 
to form each star generation: indeed, in our scenario SG and 
TG stars are created using only a relatively small fraction 
of the total mass that was used to form FG stars. In other 
words, subsequent stellar generations will be more and more 
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Figure 5. As in Figure [T] but for massive PSs. Stars of the first, second, third, and fourth generations are represented as red, 
blue, yellow, and orange circles, respectively. Panels represent: a) Formation of FG stars; b) Gas accumulation, from a mix of 
winds of FG massive stars plus pristine gas; c) First core-collapse SN explosions, which trigger the formation of SG stars; d) Gas 
accumulation, from a mix of FG SN ejecta, winds of FG and SG massive stars, and gas remaining from the formation of SG 
stars; e) Formation of TG stars, when SN explosions end; f) Gas accumulation, from winds of super- AGB and AGB stars from all 
generations; g) Formation of 4G stars; h) Present-day O-Na anticorrelation. 



affected by the evolution of lower-mass stars, whose ejecta 
will not be chemically very different from that of the original 
PS gas (panels f and g in Fig. fSJ. 
9c. Finally, a massive PS can lose the gas that has not been 
used to form stars in different ways, the most important 
mechanism probably being the interaction with larger struc- 
tures - e.g., when a PS passes through the Milky Way. 
When this happens, at first only the intracluster gas is lost 
to the GalaxyO Subsequently, however, due to the decrease 



3 According to Prie stley et alJ d201ll) . the mass evolution of the intra- 
cluster gas depends mainly on the PS mass, the PS velocity, the specific 
stellar mass loss rate, and the density of the medium where the PS is 



of the PS's potential well, the sizes of the stellar orbits in- 
crease (e.g.. iMoeckel & Bate! l20irj iDecressin et alJ 120101: 
iTrenti et alj|2010l and references therein), and in this way 
the massive PSs can more easily lose their stars as well. If 
the gas expulsion is inefficient - e.g., when the massive PS 
has a lower interaction with larger structures - the PS does 
not lose a large number of stars because the gas has been 
retained, and thus ends up as a dwarf galaxy. 



moving. However, to properly test our scenario, an extension of their 
study to larger PS masses, higher PS velocities, and denser underlying 
media would be needed. 
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Figure 6. As in Fig . El but for a massive PS after 1 3 Gyr. PGPUC 
isochrones represent the different populations presented for this 
PS: FG (red lines), SG (blue lines), TG (gold lines), and 4G (or- 
ange lines). In black lines represents a possible younger TG pop- 
ulation. Inner panel is a zoom-in at the RGB bump. 



Of course, the strength of the interaction depends on the 
mass of the host galaxy. In the case of smaller galaxies 
(e.g., the Large and Small Magellanic Clouds), the possibil- 
ity that the PS may lose mass due to interaction with the 
parent galaxy is much diminished, compared with the case 
of a galaxy like the Milky Way or M3 1 . Conversely, the op- 
posite may be true in galaxies that are more massive than 
the Milky Way, with the PSs pote ntially losing more mass 
in such massive galaxies (see also Bekki 2011). Therefore, 
clusters with similar present-day masses in the Milky Way 
and the Magellanic Clouds are likely to be associated with 
more massi ve PSs in the forme r than in the latter galax- 
ies (but see Chandar et alj|20fob . The end result is that the 
chemical inhomogeneities are expected to be less extreme 
in Magellanic Cloud clusters of a given (present-day) mass 
than in M ilky Way clusters of si milar (present-day) mass 
(see also iMucciarelli et al.l 1201 lal) - and vice-versa in the 
case of galaxies that are more massive than the Milky Way. 
This has important potential implicatio ns for the interpreta- 
tion of the integ r ated light of GCs (e . g..lKavirai et al]|2007t 
Rev et al.l 120071: iMieske et al l 12008 ^Georgiev et all 2009; 



Coelho et al.ll201 It iPeacock et al.N2.Ql lb . potentially leading 



to parent galaxy mass-related effects. 

Note that, according to our scenario, each subsequent stel- 
lar generation in co Cen need not differ substantially in age. 
This is in contrast with early claims that co Cen harbors stel- 
lar populations with an internal age spread that can reach sev- 
eral Gyr (e.g.. iHuehes & WaUersteinl l2000b iRev et al.l l2004t 
IStanford et alj2 006). It should be noted, however, that these ear- 
lier results were obtained without taking into account the possi- 
bility of different levels of He enrichment in the different popula- 
tions, which can lead to biases in the age estimates. In addition, 
it has also been recently suggested that the ages of the different 



popul ations in co Cen are indeed very similar, perhaps to within 
2 Gyr (ISollima et al.l2005b1:lPancino et al.l201 ltlD' Antona et all 
1201 ll) . That a uniform-age solution is c onsistent with t he clus - 
ter's CMD is clearly shown in Figure 8 of lSollima et al.l (l2005bb . 
where a fit to high-quality ACS@HST photometry is shown in 
the constant-age/multiple-Y scenario. 

As for the previous PS classes, in panel h of Figure we 
show schematically the expected shape of the diagnostic [Na/Fe] 
vs. [O/Fe] diagram. The corresponding HR diagrams for this 
class of PS are shown in Figure [6] 

In regard to the age issue and co Cen's CMD, one can observe 
from Figure |6]two important points: 



- The youngest SGB of the TG population with normal Y is 
not parallel to the SGB of an older FG population. Coeval 
SGBs, on the other hand, are parallel. The observations re- 
veal nearly p arallel SGBs for the d ifferent populations (see, 
e.g., Fig 5. of IVillanova et al.l l2007). thus supporting a small 
age spread in the cluster. 

- The RGB bump of the TG population shows a greater spread 
in luminosity than the one of the FG population if these pop- 
ulations have the same helium abundances. However, if the 
TG is helium-rich, the spread in luminosities are predicted 
to be similar. That the ob served spreads ar e quite similar is 
supported by Figure 10 of IRev et ail d2004 . 



As a final comment regarding the origin of the putative TG 
(i.e., the more metal-rich) component in co Cen, we note that 
it clearly corresponds to a discrete component in this cluster. 
However, it is not clear to us how several Gyr can go by af- 
ter the SN II ex plosions until such a metal-rich component is 
formed (see also [Grafton et all 1201 ll) . SN la do not appear to 
constitute a feasible solution, since the most metal-rich stars 
in co Cen a re also highly enriched in O and Na with respect 
to SG stars ( Johnso n&PilachowskfcolOl: iMarino et al.ll201 lal: 
[ Gratton et al.1 1201 ll). In fact, yields of SN II (see Fig. 8 in 
iLimongi & Ch ieffi 2003) do show an increase in Na, but not in 
O - which means that, according to our scenario, TG stars will 
have [O/Fe] similar to SG stars, but higher [Na/Fe] and [Fe/H]. 
This is precisely the behavior that is observed in Figure 19 of 
IJohnson & PilachowskH (l2010l) . 

In summary, according to our scenario, in the case of a very 
massive progenitor, we expect that SG stars are He- and Na-rich, 
but O-poor (yields from a combination between pristine gas and 
massive stars before they explode), with a small enrichment in 
Fe (contribution from a small fraction of FG SN II material). 
TG stars, on the other hand, are expected to be He-, Na-, and 
Fe-rich, but O-poor. Even though this model predicts separate 
distributions in the O-Na plane for the FG, SG and TG stars - 
a gap which becomes "bridged" when the 4G stars are formed 
(see Fig. [5] panel h) - one should bear in mind that this is only a 
schematic representation, and the real chemical evolution can be 
more complex, especially due to the fact that 4G stars are formed 
with material processed by previous generations plus a fraction 
of pristine gas (if the latter has not been completely expelled). 
In fact, as is observed in the O-Na anticorrelatio n of co Cen 
(IJohnson & Pilachowskl 120101: IMarino et all 1201 Oh . there is an 
overlap between metal-poor and metal-intermediate branches in 
the O-Na plane. However, metal-rich stars ([Fe/H] > —1.3) are 
all Na-rich, which can also be associated to some degree of He 
enrichment. 
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3.5. General Comments 

In this section, we have presented a scenario for the formation of 
multiple populations in GCs which makes several non-standard 
assumptions. The purpose of this subsection is to explain why 
these assumptions were made. 

First, we incorporate into the cluster's evolving chemistry the 
ejecta of massive stars. Most frequently, however, the contribu- 
tion of massive stars is ignored, due to the fact that their wind s 
are very fast (200 km/s < v < 2200 km/s: lLamers & Nugisl2 002). 
We decided to take these winds into account for four main rea- 
sons, namely: 

- First, massive stars are clearly good candidates to explain 
populations with high helium abundances without an in- 
crease in metals, provided their ejecta is considered before 
they explode as core-collapse SNe. This is due to the fact that 
they go through evolutionary stages durin g which their atmo- 
spher e s are comprised entirely of helium (|Limongi & Chieffil 
120071: iPastorello et al.1 l200l Uefferv & Hamannl 120101 and 
references therein). Moreover, core-collapse SNe classifica- 
tion have divided them into three groups, depending on the 
elements observed in their spectras: i) SNe Type II, whose 
spectra show H lines; ii) SNe Type lb, whose spectra shows 
He, but no H lines; and iii) SNe Type Ic, whose spectra show 
neither H nor He lines. It has been proposed that these differ- 
ences are due to the SNe progenitors' different masses, with 
the progenitors of SNe Type II being less massive than those 
of SNe Type lb, which in turn a re less massive than the SNe 
Type Ic progenitors (see, e.g., iLeonar for a recent 
review). It appears accordingly that the progenitors of SNe 
Type lb and Ic return to the ISM a large amount of helium 
before they explode. 

- In this sense, the material that arrives at the center of the PS 
is composed of a mix of pristine gas and the ejecta provided 
by the winds of massive stars. Before SNe explode, the yields 
of these winds are similar to the pristine composition early 
on, but almost completely comprised of He created via high- 
temperature H burning towards the end. Naturally, the PS 
cannot and does not select only the He-enriched wind ejecta 
to create SG stars; rather, SG stars are formed from a mix of 
pristine gas (i.e., gas that was not used to form FG stars, plus 
winds of FG stars) and He-rich winds. 

- If PSs have lost a large fraction of their mass in the course 
of their lives, and their present-day mass is of order 1 % of 
their initial mass onlyQ the escape velocity of the PS was ten 
times greater than the present-day escape velocity, thus mak- 
ing the retention of the massive stellar winds within the clus- 
ter's gravitational potential well much more feasible early 
on. Moreover, the presence of intracluster gas also plays an 
important role in decreasing the velocity of the ejected gas. 
For example, in the case of a GC - a group of stars with- 
out intracluster gas - the velocity of the ejected material de- 
creases only by the action of the GC potential well, and so, 
if the initial velocity is greater than the escape velocity, that 
ejecta will be expeled from the GC. However, in the case 
of a PS - a group of stars with intracluster gas - the veloc- 
ity of the ejecta of stars is decreased by the GC potential 
well, as well as due to the interaction with the intracluster 
gas (e.g., lDopitalll98l l : iDopita & Srriithlll986t iBrown et alJ 



ll995tlKasliwal et al.ll2005h . This means that, even if the ini- 
tial velocity of the ejecta is greater than the escape velocity, 
the ejected gas can still be retained by the PS. Naturally, this 
process would increase in importance with the mass of the 
PS, since more massive PSs have a larger amount of pristine 
gas that can interact with the ejecta. 
- If pristine gas is required to form SG stars (as implied by 
the Li and F abundances measured in GC stars; see, e.g., 
iPrantzos et al.l 12007). the latter must form before the first 
core-collapse SNe, otherwise the pristine gas has either been 
expelled or become metal-enriched. Massive star winds are 
accordingly the only viable source of gas enrichment. Even 
though one can think that SG stars will never be enriched 
in metals, it is conceivable that the amount of core-collapse 
SN explosions required to trigger the formation of these stars 
depends on the amount of gas in the center of the PS, which 
could impli es an increase of th e metallicity, as in the case of 
to Cen (e.g JPiotto et alj|2005l) . 



It should also be emphasized that stars with extreme He 
abundance may naturally have extreme abundances of other el- 
ements, such as high Na and low O. Therefore, stars at the 
high-Na, low-O end of the O-Na anticorrelation (panel h in 
Fig. [5] and panels g in Figs. [3] and Q]) are those which are 
more naturally expected to have high He ab undances. This 
is in t he same sense as recently observed by Bragagha et al. 
(l2010h . This would also be qualitatively consistent with the 
observed correlation between the presence and extent of these 
abun dance anomalies and position along t he horizontal branch 



(e.g., [C atela 



Marino et al 



& de F reitas Pachecol [1991 iGratton etal J 120101: 
201 lch . However, to more properly test this sce- 



4 iDecressin et al.1 i2007al) have pointed out that, in order to reproduce 
the high fraction of SG stars, it is indeed necessary for the GC to have 
lost close to 96% of all FG stars, unless a flat IMF for massive stars is 
adopted. 



nario, models of massive stars with mass loss for low metallic- 
ities and with time-dependent yields are required, but these are 
unfortunately not yet available in the literature. 

Regarding the ratio between FG stars and subsequent gen- 
erations of stars, the following remarks are in order. If every 
new generation is formed in the core of the PS, then the num- 
ber of FG stars decreases more than does the number of stars 
belonging to the subsequent generatio ns, due to the prefer ential 
evaporation of outer s tellar members (IDecressin et al.ll2~008l) . As 
ICarretta et al.l d2009bl) have suggested, these expelled stars can 
then become a mai n component of the ha lo field. Moreover, the 
effect presented bv lDecressin et al.l (120081) can be intensified due 
to expansion of the PS d riven by expulsion of t he gas that was 
not used to form stars dMoeckel & Batd 120101 ; IDecressin et al.l 
120101: iTrentieTalE oiO). Even though our schematic plots try to 
represent the effect of the expulsion of FG stars, these figures 
may to some extent also overestimate the number of FG stars. 

Finally, an additional case of PS could in principle exist, 
bracketing our definitions of intermediate-mass and low-mass 
PSs, where only FG and SG stars, but not TG stars, are formed. 
Such a PS would be expected to be characterized by a large gap 
in the O-Na diagram, between 0/Na-"normal" and O-poor/Na- 
rich stars. To the best of our knowledge, no such cases have been 
found in the literature. Indeed, for such a scenario to material- 
ize, a continuous mechanism responsible for removing the stel- 
lar ejecta of super-AGB and AGB stars should be present. This 
mechanism could be provided by SN la explosions; however, a 
fraction of this gas could still be retained if these explosions are 
not very close to the PS core, thus allowing the formation of TG 
stars with a mixed chemical composition. 
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4. Quantitative Estimates 

In this section we show a method to estimate the initial mass M/ 
of the PS that is expected in our scenario, on the basis of the 
present-day observed ratio between the number of SG and FG 
stars (R s p %). 

First, assuming a constant star formation efficiency e, inde- 
pendent of the chemical composition, the total mass used to form 
the FG is 



M FC =M f Xf, 

while the remaining pristine gas has a mass of 
M g = Mi X (1 - e). 

Then, the mass used to form the SG stars is 



M SG = (M g x f g 



+ Mej-fc) X e, 



(1) 



(2) 



(3) 



where f g and M ej -pG are the mass fraction of the pristine gas and 
the fraction of ejected mass of FG stars which has fallen to the 
core of the PS (and thus used up to form SG stars), respectively. 
We estimate the last term from 



M ehFC = MfQ x / x //, 



(4) 



where / is the (time-dependent) fraction of the total mass ejected 
by massive FG stars, and fi is the fraction of this ejecta that has 
arrived to the PS core - which we again assume to be entirely 
used up in the formation of SG stars. 

Now, the present-day (i.e., at time fee) number of FG stars 
is Nfc x t]fq, where Nfg is the initial number of FG stars and 
t]fg is the fraction of the Nfg stars that are still alive at time tec- 
This can be determined using the IMF (<p) as follows: 



J o.i 



(m) dm 



I (f>(m) dm, 
J o.i 



(5) 



where x is the highest mass value at time tec for the 
relevant chemical composition. For the maximium stellar 
mass, we adopted a c onservative limit for low metallicities 
dStothers & Chinll 19931). while th e lowest mass limit was chosen 
according to lKroupa et alj d!993h . 

On the other hand, the number of stars which are still bound 
to the GC after a time tec is given by 



N, 



FG—now 



J] PG X(l -S l )xNfc, 



(6) 



where S l is the fraction of FG stars that have been expelled 
from the PS by dynamical interactions with the Milky Way 
and other structures. Moreover, Nfg is related with Mpc as 
Nfg - MpG/(m) t , where (m) t is the mean stellar mass, which 
depends on the time and the chemical composition, except when 
tec = «m)o). Accordingly, 



Nfg-. 



FG-now — 



t]fg x(1 — S l) x Mpc 



(m) 

For the SG, this number is 



N : 



SG—now 



tisgX(1-S s l g )xM sc 
(m) 



(7) 



(8) 



where we have assumed that a fraction S S L G of SG stars has also 
been lost. 



Finally, the ratio between Nsc-now and Nfg-how (^fc^ wn i cn 



can be estimated from observations, is 



R 



i/s G x[(l-e)x/ g + ex/x//] 



FG 



T]FC X (1 -S L ) 



x(l 



S S L G ), 



(9) 



which interestingly does not depend directly on either M/ or 
Mgc, but is related to these quantities indirectly through S l and 
77. In eq. [9] all variables (including, naturally, R FC ^) can be es- 
timated (by models and/or by observations), except for S l and 
Sf°. Therefore, if a value for S S L G is assumedQ S l can be de- 
rived as follows: 

T] SG x\(l -e)xf g + exlxfi] 
Sl = 1 [ ^ -x(l-S s L c ), (10) 

11FG X R S F % 

Now, in order to determine the M/ of a PS once S 1 has been 
computed, we note the following: 

- Mi is related to Mpc and e through eq.[TJ 

- Mpg is related to the present-day mass of FG stars 
(M fg-,,ow), Sl and the present-day mass fraction of M F c 
(9J? FC ) by M FG-now = M FG X (1 - S X m F g- The 3R TO value 
is determined considering only stars with masses lower than 
the maximum mass (x) of stars belonging to the FG that are 
still alive, or 



J0.1 



m x (p(m) dm 



120 



m X (p(m) dm. 



(11) 



0.1 



- Mpc-now is related to Mgc and the fraction of present-day 
FG stars as M FG . now = M GC x (N F G-nowlN„ w), where N now 
is the total present-day number of stars in the GC, or equiv- 

alently, M F c-now = M cc l (l + R S FG ).E 



M, 



In summary, we find 
M cc 



m FG xex(l +R f %)x(1-Sl) 



(12) 



which relates Mj with S l, which can be estimated from the ob- 
servations using eq.[TOl(assuming a value for S s , c ). 



4.1. The Case of NGC 2808 

We now estimate S L and the M, of NGC 2808, a GC which we 
assume to be representative of intermediate-mass PSs. 



5 Observations have recently shown that ~ 2.5% of the halo stellar 
population is CN-stro ng, and thus can be asso ciated to SG stars that 
once belonged to GCs (Ma rtell & Grebell[20Toh . Moreover, the N-body 
simulations of iD'Ercole et all ([2008) show that if SG stars are formed 
in the centers of GCs (corresponding to ~ 15% of the total population), 
after 13 Gyr the fractions of FG and SG stars that are lost are ~ 99% 
and ~ 70%, respectively, with the end result that the GC population is 
dominated by SG stars (~ 80% are SG stars). 

6 This neglects the mass of the stellar remnants, which constitute ~ 
30% of the total GC mass after 12 Gyr, the exact remnant mass fraction 
depending on the detailed dynamical evolution history, metallicity, and 
the IMF (Kruijssen & Lamers 2008). On the other hand, the mass-to- 
light ratio also depends on these same para meters, varying in the rang e 
M/L v a 1 to 7 M Lg' at an age of 12 Gyr I Kruij ssen & Lamerl2 008). 
Since the uncertainty in M/Ly is much higher than that brought about 
by neglecting the stellar remnants, we decided to avoid the latter in our 
first, and admittedly rough, PS mass estimates. 
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Figure 7. Relation between the initial mass and the total mass 
loss by massive stars aft er 3 Myr. Dots are from Table 2 of 
ILimongi & Chieffil (2006), where a constant mass loss rate was 
assumed. The red line is our linear fit. 



We use a value of e — 0.4 dParmentier et alJl2008i) . while the 
values for r\ can be determined from eq. [5] if the initial mass of 
the most massive star of each generation is known; the masses of 
low-mass stars are assumed not to decrease by a large amount. 
They are calculated knowing the chemical composition, and the 
age of each star generation. Here, we assume that both genera- 
tions have the same age (tec — 12 Gyr), metallicity (Z = 0.001), 
and a-element ratio ([a/Fe] = 0.3), bu t different helium abun- 
dances (Y FG = 0.245 and Y sc = 0.345 : iPiotto etalJl2007l) . This 
implies maximum stellar masses of a 0.8 1 Mm and - 0.68 M 
for the FG and SG, respectively dValcarce et al.ll20"TTt) . Using a 
iKroupa et all (1 1993b IMF, we also find ij fg = 0.50 and i] SC = 
0.44. 

To determine /, we use the values from Table 2 of 
ILimongi & Chieff] (120071) . In this sense, as far as the present 
study is concerned, the most important conclusions are as fol- 
lows: 

- The first core-collapse SN (initial mass of 120 M ) explodes 
after 3 Myr. 

- Only stars more massive than 20 M lose a significant 
amount of mass in this short time interval. 

- We obtain a linear relation between the initial mass of the 
stars and the mass lost after 3 Myr (Ml oj ,; see Fig. |7J, which 
reads as follows (with masses given in solar units): 



M ios ,(3 Myr) = 0.88 x M initM - 18.25. 



(13) 



We assume, for simplicity, that massive stars lose their en- 
velopes at a uniform rate during their lives. 



/ can thus be estimated as 

>120 



-120 



I Mio St (m) X (f>(m) dm l I m x cf>(m) dm, (14) 
J 20 I Jo.i 



which at 3 Myr takes on a value of 0.015. In other words, after 
3 Myr stars more massive than 20 M will have ejected 1.5% 
of the total mass used to form the entire population of stars to 
which they are associated - low-mass stars included. 
Finally, eq.[10]is reduced to 



1 



5.28 x 10-' x f g + 2.64 x 10~ 3 x /, 



R 



SG 



(15) 



which depends on f g and fi, and where we have assumed, for 
simplicity, that a negligible number of SG stars has been lost. 

Due to the fact that our model predicts that SG stars should 
have low O abundances, to obtain an estimate of f g we use 
[0/Fe],„,„ and [OfFe] max , which are the minimum and maximum 
abundances of this element relative to iron that are observed in 
the present-day GC. We proceed as follows: 

- First, in a cloud of mass M u , the number of particles per unit 
mass of an element X (Nx) is 



JV X = M u x ^ 



(16) 



where Fx is the mass fraction of element X. 
Then, the ratio of the number of particles per unit mass be- 
tween O and Fe in the progenitor cloud of SG stars is 



No\ _ MgXf g x F , g + M FG x I x fi x Fq,fg 

N Fe } sc M g XfgX F F e ; g + M FG XI X fiX F Fe.FG ' 



(17) 



where M g xf g and Fx, g are the mass used and the fraction by 
mass of the element X in the pristine gas, whereas Mf G xIxfi 
and fx.FG are the mass used and the fraction by mass of the 
element X in the mass ejected by FG massive stars. 
Now, if ^Fe,g = -FVe.FG = ^Fe, which means that massive 
stars do not eject material enriched in Fe in the first 3 Myr, 
eq. [T7]can be rewritten as 



\NfJsc 



f g x(l-e)x F 0tg + ex I x fix F p jg 
\fg x (1 - e) + exl xfi\ x e x F Fe 



(18) 



where M g and M F q are related to Mj using eqs. Q] and [2] 
Relating Fo and Fp e to A^o and A^e, respectively (eq. [TBI , 
eq. [T2]is equivalent to 



1 



N Fe j sc ex(l+c) 



AWe I AW/ 



(19) 



where 



c = 



exlxfi 



Finally, eq. [19]can be transformed to relative abundances 

10 [O/Fe] SG 1 



(20) 



fX(l+c) 



^10[°/ F ^ +c xl0 [O/Fe] -], (21) 



which means that the ejected gas by FG massive stars must 
have an O abundance of 



[0/Fe] FC = log 



fx(l+c)x lQl ^ 



10 [O/Fe 



(22) 



Note that the fraction // / f g that is used to compute c is con- 
strained by 



1 -e 
exl 



10 A[O/Fe] 



- 1 



(23) 
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Figure 8. Fraction of FG stars that an intermediate-mass PS must 
lose (S l) to reproduce the present-day ratio between the number 
of SG and FG stars. Black curves represent different mass frac- 
tions of gas ejected by massive FG stars in 3 Myr (values of // in 
percentages) used to form the SG stars when e is 40% and when 
S S L G is negligible (solid lines) or 20% (dashed lines). The blue 
vertical lines represent the observed ratio between SG and FG 
stars (R s F q) esti mated from the num ber of red and blue MS stars 



(solid line, from lPiotto et al. || 2007l) and from th e O-Na anticor- 
relation (dashed line, from Carretta et al. 2009b), as observed in 
the case of NGC 2808. 



where A[0/Fe] = [0/Fe] g - [0/Fe] sc . Here we have ob- 
tained a relation between the fraction f g of primordial gas 
used to form SG stars and the fraction f] of gas ejected by 
FG stars that are used to form SG stars. Note that, by defini- 
tion, neither of these values can exceed unity. 

One can estimate the ratio in eq.l23lbv assuming that [O/Fe]^ 
and [OfFe] sc are directly given by [0/Fe] mM and [0/Fe] ml „, re- 
spectively, as indeed im plied by our model. Using the abun- 
dance ratios derived by ICarretta et al.1 (120061) in the case of 
NGC 2808 ([O/FeW = 0.4 and [0/Fe] mi „ = -1.0), we ob- 
tain // > 6180 x f g , where, as in eq. Q3] we used / = 0.015 
and e = 0.4. This implies that, if SG stars in NGC 2808 were 
formed using 100% of the ejected mass by massive FG stars 
(corresponding to 1.5% of M/), the maximum mass fraction of 
the primordial gas that can be used to form SG stars must be 
lower than 0.02%, or else the observed values of [O/Fe] cannot 
be reproducedQ 



7 Observationally, one finds that the abundances of such fragile ele- 
ments as Li and F are correlated with O (and anticorrelated with Na ; 
e.g., [P asquini et aLll2005l : ISmith et al. 2005; Denisse nkov et al J 120061 : 
Shen et al. 2010). However. IShen et al.1 conclude that the slope of the 
observed Li-O correlation is not what one would expect from simple 
pollution scenarios, claiming that the pollu ting gas mu st be somewhat 
enriched in Li. If so, and as pointed out by She n et all this would rule 
out massive stars as the main polluters, unless mechanisms can be found 
through which these massive stars can produce Li. Contrary to their con- 



Even though the derived limit for f g implies that O is com- 
pletely depleted in the gas ejected by massive FG starsj^it can be 
used to constrain the minimum value of S l depending on f g (or 
//). Then, using eq. Q3J the fraction of FG stars that have been 
expelled is 



S L > 1 -2.73 x 10" 



x///4g- 



(24) 



Figure[8]shows the minimum fraction of FG stars that NGC 2808 
has lost, according to our model. In the specific case of this GC, 
R S F G can be determined by two methods: 



- Using the MS split (Piot to et al.ll2007l) . where 63% of the 
population is in the reddest MS, while 13% is in the bluest 
MS, and both populations can be associated to FG and SG 
stars, respectively - thus implying an /?^(MS) ~ 0.20. 
However, this method is the minimum value for R SG be- 
cause, according to our scenario, an unknown fraction of 
the stars in the redder MS were formed after some tens of 
Myr from material processed by super-AGB or AGB stars 
(thus belonging to later stellar generations; see Fig. [3j. These 
young stars in the redder MS must have a chemical compo- 
sition (and in particular a helium abundance) similar to FG 
stars, but still showing abundance variations in some of the 
light elements. 

- Using the O-Na anticorrelation, where stars with low Na 
abundances belong to the FG and stars with low O abun- 
dances belong to the SG, whereas stars that belong to nei- 
ther of these groups are TG stars. Following ICarretta et al.l 
(2009b) definitions for primordial, intermediate, and extreme 
populations, R S F G can be determined dividing the fractions of 
"extreme" stars by the fraction of "primordial" stars. This 
gives an R S F G (0 - Na) = 0.36. 

In Figure [8] these R S F G ratios are represented by blue ver- 
tical lines, with the solid line corresponding to the MS CMD 
constraint, and the dashed one to the O-Na diagram constraint. 
If SG stars were formed using the total mass of massive FG 
stars (// = 100%), the minimum S L ~ 0.992 for R S F $(0 - Na), 
while for /^£(MS) the minimum S L ~ 0.986. However, if 
// = 10%, the minimum S L values are ~ 0.999 and ~ 0.998 
for R S F ^(0 - Na) and /?^(MS), respectively. Even though these 
values are high, the initial-to-final mass ratio of GCs depends on 
the orbit of each PS around the Milky Way, which can com- 
pletely destroy them if the interaction is strong enough (e.g., 
lKmiisseneTal]|2011h . 



elusions, however, the chemical evolution models bv lLind et alJ d20iTl) 
account very well for the observed Li abundances in NGC 6397, using 
massive stars between 20 and 120M o as the polluters - provided their 
ejecta are suitably diluted with pristine material. Clearly, the chemical 
evolution of Li is far from being straightforward, as the remarkab l y con- 
stant Li abundance with metallicity observed by Mo naco et al.l d2010h 
in a> Cen also shows. In addition, it is important to bear in mind that 
it is not straightforward to interpret Li abundance observations, given 
the many ill-constrained physical mechanisms that may play a rele- 
vant role in defining sur face Li abundances in present-day stars (e.g., 
iMucciarelli et aU201 lbh . 

8 In fact, massive stars do produce oxygen, albeit to a lower extent 
than they produce iron. However, the material chemically enriched in 
elements heavier than C, He, and H, in the case of the progenitors of 
SNe Type Ic, Type lb, and Type II, respectively, is not expelled via 
winds, but rather as SNe ejecta. This means that only PSs that can 
retain the ejecta of SNe can also have an enrichment in such heav- 
ier elements. This is in fact observed in the O-Na anticorrelation of 
fci Cen, when abundance data are separated into bins of metallicity 
dJohnson & Piiac howski 20Tol: lMarino et alj201 lah . 
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Figure 9. As in Figure [8] but in the case of at Cen. The blue 
vertical line represents R S F G+TG = 0.67 (see text). 



Using the range of S l values determined with // = 1 and 
eq. [12] the minimum Mj of NGC 2808 is estimated between 
3.8 and 5.9 x 1O 8 M , where we used TO .8i = 0.5040 In fact, 
the large fraction of cluster mass lost makes our scenario more 
plausible, because the escape velocity of the PS of NGC 2808 
was between ~ 860 and ~ 1060 km/s (and possibly higher, de- 
pending on how many SG stars are lost), instead of the present- 
day ~ 50 km/s dGnedin et al.ll2002h . Unfortunately, there are no 
estimates of the Mj for NGC 2808 in the literature. However, 
the orbit of this GC suggests a much more massive PS than the 
present-day GC mass, due to interaction with the Milky Way 
(ICasetti-Dinescu et al.ll2007h . 

Another important point that should be emphasized is that 
our scenario naturally explains the discrete nature of the stellar 
populations observed in NGC 2808. This is due to the fact that, 
after the first core-collapse SN explosion and the triggered for- 
mation of SG stars, there will be a period of time without star 
formation, due to subsequent FG core-collapse SNe, which are 
interspersed with SG core-collapse SNe, followed by SNe Type 
la. In fact, one should expect that this process is repeated ev- 
ery time that a new generation is born, but every time with a 
chemical composition more similar to the FG. Accordingly, the 
observational detection of the differences between subsequent 
populations (TG in this case) becomes increasingly more diffi- 
cult. 



4.2. The Case of w Cen 

As usual, the case of to Cen is more complex than for other GCs, 
due to the fact that in its CMD and in spectroscopic analyses 
alike a large numbers of populations is observed (§1). However, 



9 This assumes that SG stars have not been lost during the PS's evolu- 
tion. If, on the other hand, 70% of S G stars have been expelled from the 
PS, as suggested bv lD'Ercole et al.l ( 120081) . NGC 2808's minimum M, 
values are increased, leading to masses between 1.3 and 2.0 x 10 9 M Q . 



we used the ratios of th e different populations given in Table 4 of 
IVillanova et al.l (120071) to estimate the initial mass of this GC, as 
required to form SG and TG stars from gas previously processed 
by FG massive stars (plus a fraction of pristine gas). 

In the formation scenario for massive PSs, which is the case 
assumed for u> Cen, there is only a small difference in metallicity 
between FG and SG stars, but a sizeable difference in the helium 
abundance. As in the case of NGC 2808, we associate a large 
fraction of the reddest MS of the cluster to FG stars (57%), while 
SG stars are associated to the bluest MS (33%). However, in 
this scenario the gas ejected by massive stars that is not used to 
form SG stars is used to form TG stars - which, according to our 
scenario, will be more metal-enhanced than the SG. We associate 
these putative TG stars to the MS-a com ponent (containing 5% 
of all MS stars) of Villanova et al. (2007), which is suggested to 
be r elated to the cluster's most metal-rich RGB sequence (RGB- 
a in ISolUma et al.ll2005al) . With this in mind, we rewrite eq. [9] 
as 



R 



SC+TC 
FG 



S G-now 



+ N- 



TG-now 



N FG 



(25) 



In section I4TTI we estimated a value for / = 0.015 which 
could have been overestimated. In the case of u> Cen, we again 
use the same value, though in this case it could instead repre- 
sent an underestimate, due to the fact that after 3 Myr mas- 
sive FG stars continue to return gas to the medium that is 
used to form SG and TG stars. We also use the same value of 
tjfg estimated for NGC 2808, and assume that tjsg+tg ~ ilsc 
(nimax-SG ~ 0.68 M ), due to the small percentage of TG stars. 
Finally, using eq.|23] we obtain // > 7800 f g , which determines 
a minimum value of S l for o> Cen of 



S L > 1 



2.71 x IQ- 3 x fi / R S F G+TG . 



(26) 



In Figure [9] is shown the relation between S l and Rp ^ 70 , 
where the blue vertical line correspo n ds to R S F G+TC = 0.67, 
as obtained from the [Villanova et al.l (120071) number counts. 
In the case where // = 1, the fraction of FG stars expeled 
must be higher than 99.6% of the initial number, implying a 
Mj « 2.2 x 1O 9 M - which is between a 4 and 17 times 
higher than the estimates obtained on the basis of N-bod 



simulations (e.g., Bekki & Freemanll2003l: iMizutani et 



JN-body 

aLl l2003l: 



Wylie-de Boer & Cottrelll 120091 with favored M 7 values of 1.3, 
5.8, and 5.0 X 10 8 M o , respectively) While our estimates of M ; 
thus appear higher than the most likely (dynamical) value for this 
PS, we must stress that more detailed computations, incorporat- 
ing the ingredients described in the next subsection, are clearly 
required before we are in a position to provide more than only a 
rough estimate of a PS's Mj value. 

4.3. Required Improvements 

In previous sections we have related two GCs to a massive PS 
(u> Cen) and to an intermediate-mass PS (NGC 2808). However, 
we have not related any GC to low-mass PSs, due to the diffi- 
culty pointed out in 32.21 namely the treatment of Type la SNe. 
To solve this problem, models of GC evolution are needed which 



10 Here too we have assumed that SG stars have not been lost during 
the PS's evolution. If 70% of SG stars have been expelled from the PS 
iD'Ercole et alj |2008). w Cen's M/ value is increased to 7.4 x 10 9 M o , 
leading to overestimates in M ; by factors in the range between a; 13 and 
57 with respect to the aforementioned N-body studies. 
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must take properly into account: i) SNe explosions taking place 
at different positions inside the GCs; ii) Interactions between 
gases; and iii) The fraction of binary systems which are possi- 
ble SN Type la progenitors. 

More in general, a more precise determination of the follow- 
ing parameters is also required: 

- e: the o bserved values of this parameter vary between 0.2 
and 0.4 dParmentier et alJr2008l) . which affects estimates of 
Mi, higher values implying more mass ejected by massive 
FG stars, which can accordingly be used to form more SG 
stars. 

- IMF: As shown by (e.g.) ISkillmanl (H.008), for very low 
metallicities the IMF predicts more massive stars. However, 
a proper evaluation of IMFs for me tallicities similar to those 
observed in GCs is required (e.g.. ISollima et ai1l2007h . in- 
cluding the case of high-helium environments. 

- /: better estimates of this parameter require evolutionary 
tracks for massive stars with mass loss for low metallicities, 
along with their yields (e.g.. lHerwig|l2004l for intermediate- 
mass stars). Moreover, due to the fact that such winds are 
strongly dependent on the metallicity (low-metallicity stars 
releasing a lower amount of mass), it is possible that a trend 
between subsequent generations and the metallicity exists. 
However, it is not possible to estimate if this trend is more or 
less sharp with the metallicity, since metallicity also changes 
the other two im portant parameters used in our a pproach: the 
IMF and e (e.g.. lskillmanll2008l:lDib et~alll201 ll) . 

- Other ingredients: there are many other possible mechanims 
that may play an important role in the formation of mul- 
tiple populations with high helium abundances, but which 
we have not considere d in our analysis. Thi s includes, e.g., 
massi ve binary stars dDe Mink et all l2009h . stellar colli- 
sions (Sills & Glebbeekll201C)l) . and "localized e nrichment" 
dMarcolini et alJl2009h. amo ng others (see als o iParmentier l 
120041; ISmithl 2006. 120 10t and lde Griisl20Tol and lMTtelll201 ll 
for recent reviews). Further work incorporating these poten- 
tially important ingredients is strongly encouraged. 

Finally, a solid model for the evolution of GCs and their mul- 
tiple populations must consider the evolution of the PS structure 
after the expulsion of pristine gas, which lead s to an expansion 
of the PS and a period of subsequent star loss (Moeckel & Batd 
I2010h . Moreover, in a model where massive stars are considered 
to play an important role in the evolution of GCs, the formation 
of low-mass stars cannot be assumed as instantaneous, due to 
the fact that proto-l ow-mass stars have contract ion times of sev- 
eral tens of Myr dBernasconi & Maederlll996h . and thus could 
become c ontaminated in that period of time by the winds of ma s- 
sive stars dNewsham & Terndrupl2007HTsuiimoto et al.ll20"07h . 

5. Summary 

In this paper, we have summarized some of the most recent sce- 
narios which try to explain the formation of multiple populations 
in GCs. We have also presented a new scenario, whose chief dif- 
ference with respect to previous scenarios is the classification of 
GCs by their initial mass: i) The most massive PSs can retain 
the ejecta of massive stars (winds and core-collapse SN ejecta); 
ii) Intermediate-mass PSs can only retain the winds of massive 
stars, but not the SNe ejecta; and iii) Low-mass PSs can only 
retain the slow winds of intermediate-mass stars (see also lBekkH 
1201 lh . 

To avoid the classic argument against the possibility of using 
the ejecta of massive stars because of the high velocity of their 



winds (from hundreds to some thousands of km/s), in our sce- 
nario we postulate that: i) GCs were (much) more massive in the 
past, with escape velocities of some hundreds of km/s, a large 
percentage of FG stars and the pristine gas having been expelled 
by core-collapse SNe and/or SNe Type la, and the interaction 
with their host galaxy; and ii) After the formation of the FG of 
stars, the remaining pristine gas begins to fall to the center of 
the PS, interacting with massive stellar winds that try to escape 
the PS. An important aspect of our scenario is that it can natu- 
rally explain the discreet nature of FG and (especially) SG stellar 
populations, as have been observed in GCs. 

In general, it is exceedingly difficult to estimate the PS mass 
Mi directly from the GCs present-day mass Mac, because Mac 
depends on the GCs past dynamical history, including interac- 
tions with the Galaxy. On the other hand, our scenario allows 
order-of-magnitude estimates from different empirical estimates 
of the present-day number ratio between present-day FG and 
subsequent stellar generations that may be present in individual 
GCs. This leads to the minimum M/ values for NGC 2808 in the 
range between about 3.8 and 5.9x 10 s M Q , whereas for a> Cen we 
obtain 2.2 x 10 9 M . In the case of NGC 2808, no previous M/ 
estimate s are available, but its orbi t does suggest a massive pro- 
genitor dCasetti-Dinescu et al.ll2007h . In the case of to Cen, the 
mass implied by our scen ario is higher than obtained on the basis 
of N-body simulations dBekki & Freem an 2003; Mizutani et alj 
120031: IWvlie-de Boer & Cottrellll2009h . Further work is clearly 
required, before we are in a position to provide more than a 
rough estimate of a PS's Mj value, based on present-day observ- 
ables. 

Our scenario can also explain why the O-Na anticorrelation 
does not follow the same pattern for all GCs, where extreme 
high-Na/low-O stars are SG stars formed with material pro- 
cessed by massive stars and a fraction of the pristine gas (mas- 
sive and intermediate-mass PSs), while normal high-Na/low-O 
are SG formed with material ejected by super-AGB or AGB 
stars (low-mass PSs). Subsequent generations are formed with 
a mixture of material from pre-existing generations, where the 
last stars to form (i.e., the youngest stars) are more chemically 
similar to FG stars. 

Naturally, we are well aware that our scenario represents but 
a first "toy model" that will require extensive improvements, in- 
cluding three-dimensional hydrodynamical modeling, before it 
is in a position to provide more detailed, quantitative predictions 
regarding the nature and properties of the multiple populations 
that are now commonly observed in GCs. In this sense, detailed 
numerical inv estigations, such as the one recently carried out by 
iBekkil d201 lh . are strongly encouraged. 
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